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复 Ginzburg-Landau 方程 的 数值 模拟 
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(1- 常州 纺织 服装 职业 技术 学 院 ， 常 州 213164; 2- 南京 审计 学 院 数学 与 统计 学 院 ， 南 京 211815) 

摘 要 : 本 文 对 复 Ginzburg-Landau 方程 的 周期 边界 问题 构造 了 三 个 数值 格式 。 其 中 两 个 差分 格式 的 精 
度 分 别 为 O(7? 上 + h?)， 谱 方法 的 精度 为 O(7? + hm")， 其 中 为 方程 的 光滑 度 。 用 线性 化 分 析 的 
方法 给 出 了 格式 的 稳定 性 条 件 ， 并 给 出 了 数值 实验 。 数 值 实验 表明 ， 四 阶 紧 致 格式 的 计算 效果 最 
好 ， 既 能 达到 较 高 的 计算 精度 又 能 节省 大 量 的 计算 时 间 。 
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Ħ Ginzburg-Landau 7j f£ (CGLE) 的 动力 学 性 质 引 起 了 人 们 的 极 大 兴趣 是 。 复 Ginzburg- 
Landau 方程 在 许多 方面 有 着 重要 的 应 用 ， 如 二 维 热 对 流 ， 半 导体 激光 ， 非 线性 光纤 。 在 大 多 
数 应 用 中 ， 都 用 到 了 此 方程 的 周期 边界 条 件 。 本 文 就 研究 如 下 形式 的 CGLE 方程 

iu; 十 M *(-ciR)u?u-iu, w(at)-u(bt)  u(r,0)- uo(z), (1) 
其 中 zx € ja, b], tE (0, +00), u(z, t) 为 复 函 数 ， RÄ EXX, ^L —b-a, 则 五 为 方程 的 一 个 
周期 。 文 献 [2] 给 出 了 方程 的 一 些 动 力学 性 质 ， 并 对 于 不 同 的 初 值 条 件 画 出 了 方程 解 的 模 的 图 
像 ， 但 是 没有 给 出 具体 的 数值 格式 ， 本 文 的 主要 工作 就 是 给 出 了 计算 此 方程 的 三 个 数值 格式 。 
本 文 是 如 下 安排 的 ; 在 第 2 节 给 出 了 两 个 数值 格式 ， 第 3 节 给 出 了 格式 的 稳定 性 线性 化 分 析 ， 
第 4 节 分 析 了 一 下 方程 的 平面 波 解 ， 第 5 节 给 出 了 数值 实验 。 
2 数值 格式 
取 空 间 步 长 h = (b-a)/M, M 为 偶数 ， 时 间 步 长 为 -， 则 
Tj;=a+jh, t,—nr, j=0,1,..……,M, n=0,1,2,.…. 
本 文 所 用 的 符号 如 下 
(uj) = (一 /MT  (u)as = (1 — 2U? + v. 1) /2, 
与 文献 [1] 类似， 我 们 把 方程 (1) A tr Ilta 分 为 两 步 计算 ， 首 先 求 方程 
iu; + (1 + iR)|u[?u = 0, (2) 
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再 求解 方程 


1 
iu; 十 gurs =0. (3) 


格式 1 对 方程 (2) 精确 求解 ， 对 (3) 用 二 阶 隐 式 差分 格式 ， 得 如 下 的 分 裂 格式 


* 1 i nl2 
uš = eTl +iC+iR)|u7 *] ju, 


. ü 1 z . 
i(u?*! —uj yr t 1l Das 十 (u)z«] = 0, 小 三 0,1,2, M —1, 


此 格式 的 精度 为 O(7? + h?). 
格式 2 ”对 (3) 用 四 阶 紧 致 格式 回 ， 即 用 算 子 


C to p 


逼近 二 阶 微 商 算 子 2 BUFER 


i ; nl2 
uš = eri HOHB)Iv |] us, 


i(u*! — u2)/r 4- HU $ P pr) pigg 4uf)z-0, j20,52,.-.,M-1 
j j 4 12 7 mz j) 35% J59,1,4, , , 
此 格式 的 精度 为 O(r2 十 h^). 

格式 3 “与 文献 H] 类 似 ， 我 们 用 谱 方法 进行 计算 
T[1+i(1+iR)|u? |? 


* 一 l.u? 
uj =e uj, 


M/2-1 
up = Y e Mi T/2ge .eiiz;o-9) 4 —0,1,2,...,M — 1, 
i-— MJ/2 


其 中 


2ni 
Pi 





M-1 : 
1 : M M 
s+ *e-iu(z;-a) eo 
Za “i 2. us 25 7 88 


此 格式 的 精度 为 O(r2 +h”), m 为 方程 (1) 的 光滑 度 。 


3 ”格式 稳定 性 的 线性 化 分 析 
下 面 我 们 用 线性 化 的 分 析 方法 分 析 格 式 的 稳定 性 ， 令 好 = Ere, RARR 


uš = erii GO-HR)IE" P] , greie, 


ir ir ir n ir , iri, ir, 
ba 十 (1 + Ar e 4 = 4 T G — Shu; + gito 
其 中 r = pe BA up! = 6"+leij4， 代 入 上 式 并 化 简 得 人 +1 = G5 r6), Hp 
p 工 一 E + E cos 


Galé. r, 0) = eri RIEF] , eirlé” 2 : : 
1(6,7,0) 1-4 if — f cosó 
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格式 稳定 的 充 要 条 件 是 |G1(&,7,9)| < 1， 得 erB-Rie" 门 < 1， 即 |en|? > 去 ， 这 与 方程 自身 的 性 
质 有 关 ， 与 网 格 比 7 无关， 格式 1 绝对 稳定 。 
对 于 格式 2， 将 好 = ene? 代入 得 


u = eTl1+i(1+iR)|£”1?] .Enei78， 


(二 
ü (+ + "eam (it 3) E ( ex) 
因为 w+ = Erte, AERARII eH = G6 r0), Kp 


11. eirl£" l? $ -£ 5+ (+ it) cos0 


G3(£,r,0) = T[1- R|£"| 
Y(&r ) e Pra jo 


我 们 可 以 得 到 与 格式 1 相同 的 结论 ， 格 式 的 稳定 性 与 方程 的 自身 性 质 有 关 ， 与 网 格 比 > 无 关 。 
对 于 格式 3， 由 了 Parseval 定理 知 


M/2—-1 


lu"? = (b-a) X e - hy 


{=—M/2 


M 
= e2r[L 一 有 RE" 门 .p 》 ey = e2r[L 一 如 6 站 (b — ay (£^? 


j=1 


我 们 可 得 (e"+1)2 = e270 一 RI&" 中 (En)?， 所 以 我 们 有 与 格式 1 和 格式 2 相同 的 结论 。 


4 ”格式 的 平面 波 解 
方程 (1) 有 如 下 形式 的 平面 波 解 ， 令 





u(z, t) = Ael 于)+ot， p= s, 
k 为 整数 ， 代 入 (1) 得 
A= ig w 一 s P 
不 妨 取 4 = 项 ， 此 时 
u(z,t) = v 


其 离散 形式 为 
u(zj, t") = R 3g ) 二 ort， 
我 们 将 其 代入 格式 1 至 格式 3 验证 格式 的 精度 。 
对 于 格式 1， 我 们 有 


1 ir mos T e 
Aun (e z) ut gu = [25 uj- eere 39 ECL P 
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又 upt! ES eerun, 代入 上 式 并 作 Taylor 展开 得 








ei(e- 7 一 工 一 十 O(r2 + h2), 
所 以 有 
k? 0? 
“一 go +h), 
同 理 对 于 格式 2， 得 
kn 
w = dy Ted R* O(r? + hî). 
对 于 格式 3， 将 | 
u(z;,t") = .l QM) 
i R 
代入 得 
uj = -ági(T/Rtunr) M RS i pas 
M-1 
gl sog dem 0, s k, 
à; = 一 ro | 
Ms 2 k | M4ei(r/RHonr) e135, Do 
TS 242 
upt! ES 一 去 ei(r/R+wnr) er 证 iue 
整理 得 
OPR 
4 2p? R’ 
平面 波 解 可 以 精确 的 满足 格式 3。 
5 ”数值 实验 
下 面 我 们 用 数值 实验 来 验证 格式 的 有 效 性 。 
51 平面 波 解 


Wa = —6, b=6, k—1, h =3/8, T= 0.05, 人 = 10.0， 误 差 我 们 定义 为 


Lo = max |u} ~ u(z;,t")], 
j 


此 时 相对 误差 定义 为 L1 = Loo/A， 对 于 不 同 的 R 计 算 结 果 如 下 表 1。 

由 表 1 可 以 看 出 ， 格 式 一 的 计算 精度 是 最 差 的 ， 但 是 所 用 的 时 间 比 较 少 。 格 式 2 的 精度 是 非 
常 好 的 ， 当 环比 较 小 时 ， 甚 至 可 以 超过 格式 3 的 精度 ， 但 所 用 的 时 间 与 格式 1 相当 。 格 式 3 的 
精度 是 最 好 的 ( 当 忆 较 大 时 )， 然 而 所 用 时 间 要 比 格 式 1， 格 式 2 大 的 多 。 综 上 所 述 ， 在 要 求 一 
定 精度 的 情况 下 ， 格 式 2 是 最 好 的 ， 既 能 达到 较 高 的 精度 ， 又 能 耗费 少量 的 时 间 。 

5.2 周期 解 


R=2.0, wo(7)= R-7(10--icosz), rz €[0,2n] h-m/32, r= 0.05, 


DAAR 2 Mh lw?|， 如 图 1。 数 值 实验 5.1 和 5.2 验 证 了 我 们 的 格式 的 有 效 性 和 可 靠 性 。 
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表 1: 格式 1 至 格式 3 的 误差 比较 


Loo 

6.22(-3) 
1.00(-5) 
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1.96(-3) 
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Lı 
1.04(-3) 
1.67(-6) 
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3.27(-4) 
5.16(-7) 
2.21(-7) 
2.32(-4) 
3.65(-7) 
9.67(-8) 
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2.96(-7) 
7.96(-8) 


图 1: 格式 2 计算 的 ju?| 
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cpu-time(s) 
6.0(-3) 
1.0(-2) 
0.19 
1.0(-2) 
1.0(-2) 
0.19 
1.0(-2) 
1.0(-2) 
0.19 
1.0(-2) 
1.0(-2) 
0.18 
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Numerical Methods for the Complex Ginzburg-Landau Equation 
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Abstract: In this paper, three numerical schemes are studied for the complex Ginzburg-Landau 
equation with periodic boundary conditions. The truncation error of the two difference schemes are 
O(r? -- h*). The truncation error of the spectral method is O(7?-- ^), where m denotes the smoothness 
of the differential equation. The stability of three schemes is analyzed by a linear method, and numerical 
experiments denote that our schemes are effective and reliable. The computing results demonstrate that 
the fourth order compact scheme is the best, which can reach a high accuracy with shorter computing 
time. 
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